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I. INTRODUCTION 



The large production rate of top quarks at the Large Hadron Colhder (LHC) [1], and the 
high numbers of leptonic decays, make it a key signal for the study of precision physics in the 
Standard Model (SM). The first theoretical predictions for the unpolarised cross section at 
next-to-leading order (NLO) have been known for over twenty years [2^5} • However, recent 
years have seen substantial progress in NLO predictions for heavy quark production and 
associated processes. Spin correlations and on-shell decays of the top quarks at NLO 0, 0] 
as well as electro-weak corrections [8-14| are now well understood. One-loop amplitudes 
to higher order in the dimensional regularisation parameter, which form part of the NNLO 
prediction, have also been computed [l5[ 16|. Predictions for processes with heavy quark 



production in association with other QCD particles have also been made possible thanks to 



impressive calculations of pp tt + j 13-20], pp ttbb (2l|-24| and pp — »• tt + 2? [25[. Very 
recently, full off-shell effects of the top decays have been computed |26|, 



27 



The last few years have seen a rapid development in NLO techniques allowing the com- 
putation of new multi-leg amplitudes and cross-sections. On-shell techniques, which began 
with the development of the unitarity method |! 



29| . are quickly becoming an industry 



standard tool and, through working with gauge invariant building blocks, lead to extremely 
fast numerical evaluation for the virtual corrections to NLO processes 



30-34. Some of 



these methods have also been developed into public codes |35H37|. Together with advanced 



Feynman diagram based techniques many phenomenological predictions have been possible 



for 2 — )■ 4 processes |2l|-|27l. Last year also saw the first evaluations of 2 — )■ 5 cross- 

sections with computations of pp W + 4j ji^] and e+e~ — > 5j [45]. The extension of 
analytic unitarity techniques to deal with massive particles motivates the revisiting of the 
well known process of heavy quark production in hadron collisions. In this paper we con- 
struct compact analytic expressions for the virtual helicity amplitudes. These amplitudes 
should lead to a flexible evaluation of the NLO cross-section including both spin correla- 
tions and decays in the narrow width approximation. We demonstrate that the amplitudes 
presented here evaluate a factor of ~ 10 times faster than the analytic results of ref. j46l ] 
implemented in MCFM [47|]. Full analytic computations also offer the possibility of investi- 
gating new structures in gauge theory amplitudes. In our particular example we find new 
simplicity in the sub-leading colour contributions to the one-loop amplitudes. 

Our paper is organised as follows. We first review the decomposition of the full one-loop 
amplitude into colour ordered and primitive amplitudes which form the basic building blocks 
of our computation. We then give our notation for the spinor-helicity formalism used for 
the computation of helicity amplitudes with massive particles. In section |IV] we outline the 
unitarity and Feynman based methods employed to arrive at the compact expression. We 
give details of the pole structures and renormalisation procedure in sections IVl and IVll before 
presenting the complete set of independent helicity amplitudes in lVIII and lVIIII Some details 
of the numerical implementation are given in section IIXI before we reach our conclusions. 
An appendix listing the complete set of tree-level amplitudes is included for completeness. 



II. COLOUR ORDERING AND PRIMITIVE AMPLITUDES 



We follow the colour decomposition into primitive amplitudes as described by Bern, 
Dixon and Kosower 48|. The case of massive quarks is identical to that of the massless 



case. Although we will talk of the processes gg — )■ tt and qq — tt throughout the paper. 
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it should be noted that the amphtudes are computed with all particles considered to be 
outgoing. 



A. The gg tt Channel 

Firstly the tree-level amplitude can be written, 

4°)(l„2,3,4f) =gl J2 {T'''T^nfA\^t,2,3,4t), (1) 

P(2,3) 

while the one-loop amplitude reads, 

4')(l„2,3,4f)=y,V( Yl A^c(T»^T"3)^^4;)(l„2,3,4f) 

^ P(2,3) 

+ tr(T"^T«3)5j^4)(l,,4f;2,3)Y (2) 



where, in D = 4 — 2e dimensions, cp is defined as, 

r(i + 6)r^(i-6) 

^ (47r)2-^r(l - 2e) " ^ ' 

The one-loop amplitudes can be further decomposed into gauge invariant primitive ampli- 
tudes, 

2, 3, 4f) = APiU, 2, 3, 4f) - ^4''kl*, 2, 3, 4,-) 

- ^4klt,2,3,4,-) - ^>lfl(l„2,3,4,-), (4) 
4!)(l„4f;2,3)= 5^ |4](li,2,3,4f)+4](l„2,4f,3) + 4^J(li,2,3,4f)|, (5) 

P(2,3) ^ ^ 

where Nc is the number of colours, Nf and A'^i^ are the number of light and heavy flavours 
respectively. Explicitly the permutation set is P(2, 3) = {(2, 3), (3, 2)}. 

Performing the colour summations for the squared tree level amplitudes yields, 

2 

(6) 



j:\^f'\'-9tNc{N'c-l) J2 \Af\l„2,3At)\' - 5] Af(l„2,3,4,) 

col P(2,3) 



P(2,3) 
t\t2 _ 1 

^9tNc{N'^-i) E l4Hlt,2,3,4f)|^-^^|4°)(l„2„3„4f)|^ (7) 

P(2,3) " 

For the interference between the tree and one-loop amplitudes we find 

E4^[4°^]*=^.VA^c(A^c-1) E 4!l(lt,2,3,4,)4°)(lt,2,3,4,-)]* 

col P(2,3) 

+ (iV,2 - 1) (4;) (1,, 4f; 2, 3) - 4;;(li, 2, 3, 4,-) 

-41(1*, 3, 2,4,-)) 4°)(l„2^,3^,4f)]% (8) 
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B. The qq tt Channel 

This time the tree-level amplitude is 

(1„ 2„ 3„ ) = g', (^5^ - Af^ {U, 2„ 3„ %) , (9) 

while for the decomposition of the one-loop amplitude we find 

ih, 2„ 3„ 4f) = gtcr (iV, 5]^5^ (1,, 2„ 3„ %) + 5t^5^ ylg (1*, 2„ 3„ %)) , (10) 
with the primitives 

4; (1„ 2„ 3„ %) = 1 (1„ 2,-, 3„ 4f) - (1„ 2„ 3„ %) - ^ ' (1„ 2,-, 3„ %) 
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(Af 1 (li, 2,-, 3„ %) + ] (li, 3„ 2,-, %)) - 1 (li, 2,-, 3„ %) (11) 



and 



(It, 2,-, 3„ %) = 1 (It, 3„ 2,-, %) - (It, 3g, 2,-, %) - ^4^1 (It, 3„ 2,-, %) 

+ -i^ (4^^ (It, 3„ 2,-, 4t-) + 4'' (It, 2,-, 3„ 4t)) + -^4'^] (It, 3„ 2„ 4t) . (12) 

c c 

Evaluating the colour ordered squared tree-level amplitude yields then 

E = ^'(^' - 1) l4°Hl*,2,-,3„4t-)P , (13) 

col 

where as the interference term becomes 

^4^)40)]* = g^'^crN^iN! - 1) 4J(lt, 2„ 3„ 4t-) 4°)(lt, 2„ 3„ 4,)]* . (14) 

col 

III. SPINOR-HELICITY FORMALISM 

For massless particles it is possible to completely decompose all momenta into a basis of 
two component Weyl spinors since, 

p^ = l{p\Y\p]. (15) 
The polarisation vectors and fermion wave functions then fit easily into a helicity basis: 

u+{p) = \p), M-(p) = b], (16) 

The situation for massive momenta is a little more complicated since f[T^ is no longer valid 
and the definition of unique helicity state is no longer possible as one can always find a 



boost, through the rest frame, such that the hehcity state is flipped. However, given a 
massive momentum P, one can construct a hehcity basis by defining a specific reference 
frame with respect to an arbitrary massless vector rj (49| . 



where a (3 



= aP^ + /3r]p 

2 

{P^'np)VipP^\ ' "^^^ ^ ^^"^ ^ spinors can then be defined by: 



M±(P,m; P\7^p) 
v±{P,m;P\r]p) 



{^ + m)\7]p±) 
{^-m)\r]p±) 

y^{P'T\VP±)' 



u±{P,m-P\i^p) 
v±{P,m;P\r]p) 



{VPT\{^ + m) 

,/^{r]pT\P^±)' 
{riPT\i^-m) 



^{7]pT\P'±)' 



(18) 

(19) 
(20) 



in the following we will set a = 1 , /3 
although a choice of. 



which corresponds to the choice of reference 5^ 



a 



1 + 



4m^ 



/3 



4m^ 



VP 



(21) 



corresponds to the basis used in reference 51, H^] which also has a convenient massless 
limit. A specific choice of the reference vector that allows us to match the four component 
representation of ref. 19| is presented in Appendix |Bl Keeping the reference vectors r] 
arbitrary allows us to relate the heavy quark helicity states: 



U-{P,m; P\t]p) 
f_(P, m; rjp) 



m 



u+{P,m;r]p,P^ 



{P^rjp) 



m 



v+{P,m;r]p,P^ 



(22) 
(23) 



The spinor-helicity formalism is well suited for numerical evaluation and had been imple- 
mented for mathematica within the SOM package [53| . 



IV. METHODS 

Each of our primitive one-loop amplitudes can be written using the standard basis of 
scalar integral functions and a rational term up to higher order terms in the dimensional 
regularisation parameter e, 

/I4 [L, Z, O, ^) — <-^4.i|2|3|4 -'4;1|2|3|4 
2 3 i-1 

. V V V rt^i /-[^i 

^ ^3;i...j-l|j...fc-l|fc...i-l ^•i\i...i-\\j...k-\\k...i-\ 

+ E E c'S.,-!,...-! 4£,-,,...- + ci-i/ri + fii-i + (24) 

1=1 j=i+l 

The indices k, 1} denote the momenta at the external legs of the integrals and the sums 
run over all cyclic partitions. The dependence of the coefficients and the rational term on 
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the helicity states of the external particles has been suppressed above. Figures show 
the explicit basis integrals for the five independent primitive amplitudes. Heavy and light 
flavour fermion loop contributions have also been included. 

The computation was performed in two distinct parts. The first used generalised unitarity 
to compute the compact expressions for the coefficients of the scalar integrals. A Feynman 
diagram based approach was then taken to obtain compact forms for the remaining tadpoles 
and rational terms. 

The computation was performed in the Four-Dimensional-Helicity (FDH) scheme and 
mass renormalisation was performed using the on-shell scheme. 



A. Generalised Unitarity 



The use of the generalised unitarity |32|, |5J-|58| , extended to massive propag ators 
GjJ, is the primary reason we are able to obtain compact representations for the one- loop 
helicity amplitudes. The analytic techniques we have employed are by now well covered in 
the existing literature and have been used extensively in recent analytic computations of 
pp ^ H + 2j [621-64 1 and pp — >■ Wb b (65| production. We refer the reader to recent reviews 
on the topic for further details 66|, l67|. However, since analytic computations for massive 
amplitudes are covered to a lesser extent, we outline some of the techniques specific to our 
process. 
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FIG. 1. The 7 cut diagrams contributing to the left-moving primitive amplitude in gg 
Red dotted lines represent massive fermions, plain lines represent gluons. 
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FIG. 2. The 7 cut diagrams contributing to the right-moving primitive amplitude in gg 
Red dotted lines represent massive fermions, plain lines represent gluons. 
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D- dimensional cutting procedures |32|,|59|,l68|,|69| could also be applied to the computation 
of rational terms in massive amplitudes as described in ref. [l9|. However, as we discuss 
briefly in section IIV A3t the procedure can lead to large intermediate expressions when 
followed analytically. On-shell recursion relations offer an attractive alternative to obtain 
compact expressions directly (70l-[73j. At the present time such techniques have yet to be 
extended to the massive case. A Feynman diagram approach therefore offers a simple way 
to obtain compact analytic expressions for our process. 
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FIG. 3. The 7 cut diagrams contributing to the sub-leading colour primitive amplitude in gg 
Red dotted lines represent massive fermions, plain lines represent gluons. amplitude. 
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FIG. 4. The 7 cut diagrams contributing to the leading colour primitive amplitude in qq — )• tt. 
Red dotted lines represent massive fermions, plain lines represent gluons and blue dotted lines 
represent massless fermions. 
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FIG. 5. The 4 cut diagrams contributing to the sub-leading colour primitive amplitude in qq — )• tt. 
Red dotted lines represent massive fermions, plain lines represent gluons and blue dotted lines 
represent massless fermions. 
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1. Three-Mass Triangle Coefficients 



In this section we give an explicit example of the computation for the three-(external)- 
mass triangle coefficient C3;i|23|4. This gives a good example of a computation that is specific 

to the massive case. We consider the + H h helicity configuration in the gg — > ti since 

the + + ++ configuration is zero. Following Forde's method 57|, the triple cut can be 



parametrised by a one- dimensional complex contour integral over t. The on-shell constraints 
for the loop momentum are solved in general leaving the integrand as a rational function 
of t. The scalar triangle coefficient is then given as the boundary value of this triple cut 
integrand: 



C^S|l|23(lt^2^3-,4±) = -i5^inf,[C'[fJ|,|,3(l+ 2+,3-,4±)] 



2 

7± 



(25) 



to 



h\=± h2=± /i3=± 

where inft is computed by taking a Taylor expansion around t = oo. The complex parameter 
t appears in the parametrisation of the (on-shell) loop momentum. We construct a spinor 
basis for the complex loop momentum using two massless vectors and In our case 
the massless vectors can be constructed from p4 and pi such that, 

Kl- = ^i^P'-^y') ^ (27) 
kY = ^^^^r"^>'\ (28) 



where 7± = pi ■ p4 ± \/ {pi ■ Pa)"^ — rn^- The on-shell loop momentum li can then be written: 



K = a {kI''' - Kl''^) + '-{Kll^^lK',] - ^^^^{K',\r\K\], (29) 
where 



2 

m 



7 — 



(30) 



The procedure from this point is rather straightforward. For illustrative purposes we find 
it convenient to expand the helicity sum in eq. fl2^ and explicitly remove dependence of 
the internal reference vector of the massive spinor. This results in considerably shorter 
expressions for the integrand though in general we find that specific choices of this vector 
yield comparable sized expressions after expansion. Inserting the relevant tree amplitudes 
from Appendix |X] we find (the reference vectors for the internal gluons were chosen as 

^i", = = m\r\K\]): 

_,-/=<W /1+ 9+ q- A + \ _ 
*'-^3;4|l|23'^-^t ' ' '^t J ~ 

m(X^|l|/2] {{v4Vi){hm\] + {hVi){hV4)[Klh]) ihS)^ 



(7744^) {K',h) ikh) {k2) {k3) (l^r/i) (23) [K\l 



3j 



m(/33)^(ir^|4|/3] ((r/4r/i)(/2|l|i^t] + {l2Vi){hv,)[K\h]) 
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We are then left to feed in the parametrisation of fl29p . expand around t = oo and extract 
the coefficient of t^. This analytic form still contains the dependence on 7-1- and we must still 
perform the sum before we arrive at the full scalar triangle coefficient. In principle this can 
be done numerically but in this case we find the analytic form after the t-expansion is rather 
lengthy and we can improve the situation considerably by performing the sum analytically. 
The procedure is straightforward but tends to lead us through relatively large intermediate 
expressions. We begin by removing all "flatted" spinors in place of an explicit polynomial 
in 7. To do this, we make sure not to introduce spurious denominators by making use of: 

(2|/^J|3] = -^^(2|1|3]. (32) 

The sum of the two solutions 7± is obtained from this polynomial form by rearrangement 
into a function of: 

7+ + 7- = 2pi ■ Pa, 7+7- = (33) 

We are then left to partial fraction the spinor products. After the dust clears, we are left 
with a relatively compact form: 

^'-"3;4|l|23l-^t ' ^ ' ^^t ) ~ 

6((r/4r7i)(23) - 2(2r/4)(377i))(3|l|2]W 
/3^S232(r7il^)(r/44^)(23) 
2((2r]4)(3r/i) - {r]Ar]i){'2?>)) (^23 + 2(2|1|2]) (3|l|2]m3 

/32s23(r7il^)(r/44^)(23)(2|l|3] 
2 (2m^ + (2|1|2]) (2r/i)(2774)(3|l|2]m=^ (523 + 2(2|1|2]) (2r7i)(3r/4)[32]m=^ 
/32(ryil^)(r744^)(23)(2|l|3]2 + {r,,l^) {r^,A^) {2\l\?>f 

(2r7i)(2r/4)(3|l|2]^m (3r/i)(3r/4)(3|l|2]m 312(7747/1) (3 |l|2]m 
2(r7il^)(r744^)(23)(2|l|3] 2(771!^) (r/44^) (23) (r/il^)(r744^)(2|l|3] 
((3r/i)(3r74)(2|l|3] - (2r/i)(2r/4)(3|l|2])(3|l|2]m 

4/32(r/il^)(r744^)(23)(2|l|3] 
3(8(2|l|2K + g^3)((3r/i)(3r74)(2|l|3]-(2r7i)(2r/4)(3|l|2])(3|l|2]m 

4/3%232(ml^)(r/44^)(23)(2|l|3] 
(2r7i)(2r/4)(2|l|2](3|l|2][32]m (2m^ + (3|1|3]) (3r/i)(3r74)[32]m 

2(r/il^)(r744^)(2|l|3]2 + 2(77il^)(r/44^)(2|l|3] 
(2m^ + (2 1 1 1 2] ) 323 (2r7i ) (2r/4) (2 1 1 1 2] [32]m _ l3^s,2sUTiArii)m 

2(r7il^)(r744^)(2|l|3]3 2(r7il^)(r744^)(2|l|3]2' ^ > 



+ 



+ 



where 13 = \Jl-^- Although one may not consider this form particularly elegant we 
note that it is much shorter than expression before the 7 sum. It is also among the most 
complicated coefficients that we encountered in the computation. 

2. Bubble Coefficients 

In this section we apply the Taylor expansion method to the computation of the bubble 
coefficient C^^\^^^J^lf ,2'^ ^'i'^ ^A'^). The on-shell constraints can be solved in general leaving 
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two free complex parameters which we label t and y. Owing to the two dimensional complex 
integration the pole structure is rather more involved than the triple cut considered above 
and we must consider triple cut contributions as well as the double cut in order compute 
the full coefficient 0. Our coefficient can therefore be written as: 

C2;12;34(1^2+,3+,4t) = -zinfi[inf,[C'2;i2|34(1^2+,3+,4)]] 

- ^Ei^f*[^3;i2|3|4(l,^2+,3+,4)]] - i5^inf,[C73;i|2|34(1^2+,3+,4±)]], (35) 
y± y± 

where we choose to parametrise the loop momentum as, 

with K\ = Pi — i^^Ps- In fact, with this choice we find the second term above vanishes. 
This is particularly convenient analytically, however, for a direct numerical evaluation, an 
independent massless vector would yield more stable results. Using the compact expressions 
for the tree level amplitudes and summing over the internal helicities gives the integrands 
for the two non-zero contributions to be, 

C'2;12|34(1^2+,3+,4±) = 

(2|l|2]2(13)[/i2][/i3](2/i)(r/il^)(r744^) ^ (2|l|2](/i2)(/i3)(r7il^)(r/44^) ' ^ ^ 

and, 

2m[/i2]3(/i|4|3](3|l/3|3)([/i3](/ir/i)(/ir/4) + {h\A\3]{r]^r]^)) 



C3;l|2|34(lt ,2 ,3 , 4^- 



(2|l|2][/i3](/i3)(3|/3|/i](r7il^)M^) 



tm^3\l\2f[h3]{h2y{r]ir]^) ^ ^m3[/i3]((77i|/3|2](r/4|/3|2] - [2|/3l|2](ryir/4)) 



(2|l|2](/i3)(2|/34|3)(3|/34|3)(r/il^)(r/44^) (2|l|2][/i2](/i3)(3|/3|/i](r/il^)(r/44^) ' 

(38) 

Just as in the three-mass triangle we find a considerable benefit from taking the extra effort 
to find a closed form for the triangle contributions which are free of square roots. Many 
terms cancel between the double and triple cuts so that the final result is simply: 



.[L] ^ m=^[32]((7^i|(l + 2)(2 + 3)|r/4))- W 

^2;i2|34l^t ' ' *^ (2|l|2]2(r/il^)(r?44^)(23) 



3. Tadpole and On-Shell Bubble Coefficients 



The computation of the tadpole coefficients directly from unitarity is complicated by the 
fact that the wave-function renormalisation contributions cause the double cuts to diverge. 
A way around this problem has been introduced in the context of a numerical application 



of /^-dimensional generalised unitarity in reference [19||. The method should apply equally 



^ The alternative method of spinor integration produces an identical set of poles as can be seen from 



Mastrolia's evaluation of the double cut via Stoke's theorem 



58| 
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well to analytic evaluation, but since it breaks gauge invariance, it can lead to large inter- 
mediate expressions. Such cuts also require the six-point tree-level amplitudes rendering the 
computation more difficult. Nevertheless, methods using spinor integration technique have 



been proposed [75|, 126 



Both the tadpole and on-shell bubble coefficients combine to give the coefficient of the 
log(m^) contribution to the full amplitude. The coefficient of this logarithm is completely 
fixed by the universal IR constraints once combined with the knowledge of the cuts consid- 
ered in the previous section. With this in mind, we rearrange the integral basis as, 

n o A\ _ ry[X] t[X] 
/I4 {1, Z, O, ^) — '-^4.i|2|3|4 -'4;1|2|3|4 

2 3 i-1 

, V V V rt^i /"[^i 

^ Z-^ Z\i...j-\\j...k-\\k...%-\ ^?,;i...j-\\j...k-\\k...i-\ 

i=l j=j+l A:=j+1 

+ E E 4J.'..,-,„...... H- + + (40) 

i=\ j=i+2 

For the gluon fusion channel the universal poles structure implies, 

(41) 

(42) 
(43) 

(44) 
(45) 
(46) 

In this case we returned to a Feynman based computation of the rational terms and remaining 
tadpole contributions. Owing to the simplified form of the reduction algorithm, and the fact 
that such contributions are independently gauge invariant, the approach gives a simple way 
to reach a compact form for the full amplitude. The final results for the tadpole coefficients 
turn out to be remarkably simple and suggestive that an approach based on matching with 
the universal IR and UV structure, as proposed in ref. (sof, would generalise to all massive 
amplitudes. A complete description, however, remains for future study. 





'-'2;12|34 


'-'2;23|41 


+ 


2 ^ 


r\-^\ _ 


^2;12|34 


^2;23|41 


+ 


2 ^ 


r^H] _ 

^2;m? ~ 


^2;23|41 









and the quark annihilation channel. 



^2;'m? ~ 


'-'2;12|34 


^2;12|34 


+ lAf 


^2;m? ~ 


~^2;23|41 


^[slcfi] 
^2;23|41 


/I4 , 


r^H] _ 

^2;m? ~ 


^2;23|41 


3 ^ 





B. Feynman Diagram based approach 

To compute rational contribution and provide numerical cross-checks for cut-constructible 
parts of the amplitude we have performed two independent calculations based on the tradi- 
tional Feynman diagram approach. 

The diagrams were generated with DIANA [t^] and then further processed analytically 
with two independent FORM j?^ codes, to generate tensor integral representations of the 
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colour ordered one-loop amplitude. The tensor integrals were reduced to scalar boxes, tri- 
angles and bubbles in 4 — 2e dimensions using standard Passarino-Veltman reduction in 
one case and to scalar integrals in shifted dimension according to [79i-t82i] in the second 
case. We performed calculations using both Dirac-spinor and spinor helicity methods in 
Four-Dimensional-Helicity (FDH) and 't Hooft-Veltman (HV) schemes, which allowed us to 
verify that scheme dependence is in agreement with predictions of 8^ . 



To maintain amplitude invariance under gauge transformations of gluon fields it is nec- 
essary to perform finite renormalisation of the heavy quark mass ^B]- This is done by 
including mass counter-term diagrams with the mass renormalisation defined by, 

m = (1 + 5Z„)m^, (47) 



where 

<5Z^ = -glcrCF ( - + Slog ( ^ ) + 5 ) . (48) 



3 , ...Jl^l 



V. POLE STRUCTURE 

We have verified that our amplitudes satisfy the well known universal Infra-Red and 



Ultra- Violet pole structures 



83[. These can be broken down into the contributions from 
each primitive amplitude [l9 . Representing the divergent parts of the amplitude by the 
function V^^^ we can write, 

Afl = U[^lAf + (49) 
For the gluon channel we find (/3 as defined after eq. 



..,,.,3,4,-,^-lH-l-i,o.(||L)-l,o.(-^), (50) 

VI Xlt,2,3,4t-) = — log — — . (51) 



For the sub-leading colour we list the poles of the full colour ordered amplitude which is 



proportional to Af (1^, 2, 3, 4f) + {U, 3, 2, 4f) = A^ {U, 2^, 3^, 4? 



Ut-l(l. 2 3 4,) =i«log + i^log 

U„A^,4,J ^^^^ ^°S^(2|l|2pJ + S2,e H(3|l|3]2 

The poles of full colour and helicity summed interference with the tree level can be written 
in terms of spin correlated Born amplitudes. 



c,h c,h 



2 J2 ^^'^ ■ i^^'^V =2iVc 3)|2 2, 3, 4f) + |^(°)(3, 2)1^ l^W(l,, 3, 2, 4f) 

c,h ^ 

^i4°¥^''^i-i4';irv^f^'^j, (53) 
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where 



j:\Af\2,3)\' = 



c,h 



h 



and Ec.?«l^?H2,3)p + \Af\3,2)\'^ = '£c,h\^f^\'^- Turning our attention to the quark 
channel the analogous structure can be written as: 



+ /3 



(55) 

(56) 



\/['"lfl. 2-3 4.-) = - + — - - loe 

V^l^'^] li,2,-,3„4f = log ^ log T 

e e V ■523/ e .523P 

V^[^l(l„2,-,3„4,-) = V-[/l(l„2,-,3„4,-) = |, (57) 
2Y,A^l^-[ATr = 2Y,\AT?lNM"'\lu2,,3,At) - {Nf + NhW^^^ 

c,h c,h 

- |- (V['^](li,2„3„4f) - r['^l(l„3,-,2„4,-)) - i-V[^'^](l„2„3„4f)|. (58) 



VI. RENORMALISATION AND SCHEME DEPENDENCE 

The one-loop amplitudes presented in the previous sections still contain UV divergences 
which need to be renormalised. At this point we want to remind the reader that we already 
included the mass renormalisation for the top quark mass, which we define using the pole 
scheme. Thus we are left with the wave function renormalisation and the running of the 
strong coupling constant. For the former we are using an on-shell prescription which takes 
into account all self energy contributions, whereas the coupling constant is renormalised 
in the MS scheme. To decouple the top quark from the running of a^, one subtracts the 
diagrams with a top quark in the loop at zero momentum transfer. Following the 8^, |85| we 
define the renormalisation constants for the gluon and fermion fields and that of the strong 
coupling. 



G« , = (1 + 6Zal2)G,,, , = (1 + 5^g/2)^Q , {g^ = (1 + 5Z,^)gl , (59) 

where, in the FDH scheme. 



5Zq = 5Z^, 5^c = -^f4) ' (60) 

3e \m^' 



13 



with Cp = Details of the scheme dependence of the couphng constant renormahsation 

can be found in reference 86|. Summing up all contributions the renormalised amplitudes 
are given by, 

4'^'""^" {lu 2, 3, 4f) = (li, 2, 3, 4f) + {5Zg + 5Zq + 5Zg^) Af (1*, 2, 3, 4f) 
= 4^ (li, 2, 3, 4f) - (7.'cr4°^ (1*, 2, 3, 4f) x 

J ^-^ + C^(- + 31og(£f +5 , (62) 



3e 3 \e 

(l),ren/-. ^ „ ,1 _\ _ ^(l)/'! O Q /I -^ i ,'X'7_ i X '7 ^ /I W 



1 

At^'^"^ {lu 2,-, 3g, 4f) = (It, 2,-, 3„ 4f) + (5Zq + 5Z,J AT (1*, 2,-, 3„ 4^ 
= 4^ (It, 2,-, 3„ 4f) - ^2cp4°) (1„ 2,-, 3„ 4f) x 

riliV,-2iV; AT, , ^ /3 , , ,\ 2 2,_//i| 



V 



+ - + 31og^ +5 -f-^log^ . (63) 



3e 3 ' V e ° V m? J / 3e 3 \ rn? 



To convert these renormalised amplitudes to the 't Hooft-Veltman (HV) scheme we follow 
the well known universal structure (83l-[86l|. For our amplitudes this can be summarised as, 

4^)' [^"'l (It, 2, 3, 4t) = 4') (It, 2, 3, 4t) - glcrCpAf (It, 2, 3, 4t) , (64) 



AT' (li, 2,-, 3„ 4t) = 4'^ (It, 2,-, 3„ 4t) - gtcv [2Cf - y J 4 ' (1^, 2,-, 3„ 4?) . (65) 
After renormahsation, it becomes, 

^(1). [HV],en (^^^ 2, 3, 4t-) = 4'^"^" (1*, 2, 3, 4t-) - glcr^AT (It, 2, 3, 4,) , (66) 

^(1). [HV],re„ (^^^ 3,, 4t-) = 4'^''''" (1*, 2„ 3„ 4t-) - glcrCpAf (It, 2„ 3„ 4,) . (67) 



VII. TREE-LEVEL AMPLITUDES 



The tree-level amplitudes can be computed using BCFW recursion relations |87l . 
together with the spinor conventions of the previous section. 

The three-point amplitudes are calculated directly from the Feynman vertex using arbi- 
trary reference vectors. As an example we consider the 



/l3(l+, 2^^ , 3^) = ^ MT(1, m; (2, 0^+(3, m; 3^ r/3), (68) 



V2 

which can be expanded to give 



^ (r^i|(/ + m)^,^ (2,0(0 -m)|r^3) 



^3(1+ 2^^ 3,+) = l\ i [ (69) 

The tree-level helicity amplitudes, expressed in terms of spinor products, are listed in Ap- 
pendix |A1 



VIII. ONE-LOOP AMPLITUDES 

In this section we present a complete set of one-loop helicity amplitudes needed for tt 
production at hadron colliders. 
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A. Notation and Integral Functions 

We define the general scalar integral as H 



/^l^r(i - 26) r J, 

z(7r)-T(l - e)2r(l + e)J Utii^ ' - mj 



{ki} are the sums of external momenta entering each of the n vertices of the graph. We will 
denote the basis integrals for our helicity amplitudes as, 



-^4;l|2|3|4 


= h{m^, 0, 0, m^, si2, S23, m^, 0, 0, 0) 


(71) 


jni 

4;1|2|3|4 


= hirn^, 0, 0, m^, Si2, S23, 0, m^, m^, m^) 


(72) 


-^4;1|2|4|3 


= h{m^, 0, m^, 0, S12, Si3, 0, m^, m^, 0) 


(73) 


-^3;12|3|4 


= /3(sl2,0,m^m^0,0) 


(74) 


-'3;12|3|4 


= /3(si2,0,m^0,m^m^) 


(75) 


-^3;13|2|4 




(76) 


-'3;13|2|4 


= /3(si3,0,m^0,m^m2) 


(77) 


-^3;2|3|41 


= J3(S23, 0,0, 0,0,0) 


(78) 


-'3;2|3|41 


= /3(s23,0,0,m^m^m2) 


(79) 


-^3;1|23|4 


= /3(s23,^«^"^^0,0,m^) 


(80) 


-'3;1|23|4 


= /3(s23, 0). 


(81) 



As described in section IIV A3I we find it convenient to move all log(m^) dependence into 
the on-shell bubble. Therefore our helicity amplitudes are written in terms of the following 
two-point functions. 



2;12 



/2(si2,0,m2)-J2(m2,0,m2), (82) 

= I2{s23,m\m^) - l2{0,m\m^) , (83) 

/2;23 = /2(S23, 0, 0) - him\ 0, m^) + 2 + 0(e) , (84) 

l2-m = l2{m\0,m^), (85) 

l^^ = h{0,m\m^), (86) 

and the tadpole function is removed via 

I^{m^) = [him^, 0, m^) - l) + (9(e) . (87) 

Finite box functions are used to make the IR poles explicit: 



-^4;1|2|3|4 — -'4;1|2|3|4 ~ (2|1|2] 



^ This follows the conventions of qcdloop which was used for numerical evaluations [8 
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We also find benefits in cancellation of spurious poles that can be made explicit through the 
use of higher dimensional integral functions: 

^Sl'W ^ 2(2|1|3](3|1|2] (-(2|l|2]W4;i|2|3|4 

+ ((2|1|2] + 2m2) S23/3;l|23|4 + 2(2 1 1 |2]2/3.i2|3|4 + (2| 1 |2]s23/3;2|3|4l) , (89) 
^^"^1121314 = 2(2|1|3](3|1|2] (-^2|l|2]2s23^2jm^2|3|4 

+ (2|l|2]s23/3%7i|23|4 + 2(2|1|2] ((2|1|2] + 2m') I^,,^,^, + S2s ((2|1|2] + 2m') I^,^,^,,). (90) 

Since all the expressions are quoted for the same configuration of heavy quark helicities 
we find it convenient to factor out an overall normalisation (r7il'')(7^44'') and define, 

if = -i(77il^)(7744Vf . (91) 

B. Primitive amplitudes for gg tt 

-t{va'){V44.')A['^ (1+ 2+,3+,4) = 

/4;l|2|3|4 (m (771774) [23] ) - F2;i2(^ MWW J 

1 - (92) 

+ ^ir(l+,2+,3+,4)x(/2^„-l) 

m[23]((77i7;4)(2|l|2] + (?7i|i^i2i^23|774)) m((77i774)(2|l|2] - (27/i)(3774)[23]) \ 
2(23)(2|1|2] 3(23)2 J 

i{Vil'){V4^')AP{ltXXAt)^ 
_ / (377i)(3774)[23]W(2m2 + (2|l|2]) _ m3(277i)(2774)[23]2 (2m2 + (2|1|2]) 

^4;i|2|3|4^ 2(23)(3|1|2] 2(23)(2|1|3] 

m^[23] ((?7i?74) {2m' - S23) + (?7i|^i2i^23|?74)) 
(23) 

_ ,™ f {2m'-S2s) {vm)[23]m^ 



'3;1|23|4 



(23)(2|1|2] 



, m3[23](2(77ir^4)(23)+4(27^4)(3r;i)) , 7n3[23](277i)(2774)(2|l|2] 

+ -'s-i9iai4 n r 



+ 



(23)2 ^ (23)2(2|1|3] 

m^[23] ((37yi)(3774)(2|l|3] - (27yi)(27;4)(3|l|2]) Tn^ [23] (3771) (3774) (2|1|2] \ 
(23)2(2|1|2] (23)2(3|1|2] J 



(93) 



_ / 7n3(277i)(2774)[23]2 _ (3771) (3774) [23] 2 77^-^(771774) [23] \ 

3;2|3|4i\^ 2(23)(2|1|3] 2(23)(3|1|2] (23) J 



2:12 



7r73[23] (2,Si2(7/l7/i) - (//i|A'l2A'23|//4)) 



-lif (l^2+,3^4±) (/2;„.-l)-( 
where Kij ^ Pi + pj. 



(23)(2|1|2]2 

m[23] ((77i|Xi2ir23|r/4) + (771774) (2|1| 2] ) \ 
2(23)(2|1|2] ; 
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if (l,^2^3-,4±) X - ..3(2|l|2]74;l|2|3|4 + ^h^m - ^^ji"^]^'"' 

/ 2 f 4 \ / 3m(3|l|2] ((3r?i)(3^4)(2|l|3] - (2r?l)(2r?4)(3|l|2]) ^^ 

+ (^/3;l|23|4 + -/2;23 - "J X 4(23) (2| 1 13]/3^ J 

/ 4 \ ^6m=^(2|l|2](3|l|2] ((377i)(3774>(2|l|3] - (277i)(2,74>(3|l|2]) 

+ I -'3;l|23|4 1 X 



S2sJ " V (23)(2|l|3]si3^4 
6m3(3|l|2]2((2r74)(3r7i) + (2771) (37^4)) \ m(3r7i)(3774)(3|l|2][23] ^ 

^3)5^^ ; 2 ^^^^i^i^i^ 

, m(3r7i)(3^4)(3|l|2][23] ^ m(3r/i)(3^4)(3|l|2] ^ n 

+ 2(2^"^^ (23) ^'^^'I'l' + l^'^'' ~ ^) 

m(2r;i)(3r/4)(3|l|2](2(2|l|2] + g23) 6-2. / m3(3r/i)(3774)(3|l|2][23] 

(23)(2|l|3]s23/32 ^ (2|1|2]2 

m(277i)(2774)(3|l|2]((2|l|2]s23+(2|l|3](3|l|2]) 2m^mr],) {3\1\2] 

(23)(2|1|3]2 (2|1|3] 
m(r?i7?4)(3|l|2]((2|l|3](3|l|2]-(2|l|2]^) 2m3(2r?i)(3r/4)(3|l|2][23] \ 
(2|1|2](2|1|3] ^ (2|1|2](2|1|3] J 

/ m(3|l|2]^ ((^i774)(23) + 2(2771) (3774)) m(277i)(3774)(3|l|2]((2|l|2] + 2m') 
'''l''I^V (2|1|2](23) (23)(2|1|3]^2 

m(3?7i)(3?74)(3|l|2]((2|l|2]^^(2|l|2]g23 + 2m^523) 7^^(771774) (3 11|2] 
2(23)(2|1|2]2 (2|1|3] 

7n(277i)(37;4)(3|l|2]((2|l|2] + 2523) m3(27;i)(27/4)(3|l|2][23](2(2|l|2] + 523) 

(23)(2|1|3] + (2|l|3]2s23/32 

7n(277i)(2774)(3|l|2]^ 7n(3|l|2][23] ((27;i)(2774)(3|l|2] - (37;i)(37;4)(2|l|3]) \ 
2(23)(2|1|3] 4(2|l|3]s23/32 J 

7n(3|l|2] ((771774) (2|1|2] - (2771) (3774) [23]) 7ri(27/i)(27/4)(3|l|2] ((2|1|2] + 2771^) 



2:23 



(2|1|3]S23 (23)(2|l|3]2/32 

12m^ + 523/3^ 777(3|l|2]^((277i)(3774) + (377i)(2774)) 

(23)si3^2 

12771^ + 523^^ 777(2|1|2](3|1|2] ((277i) (2774) (3| 1|2] - (377i)(3774)(2|l|3]) \ 

S23/32 (23)(2|l|3]5i3/32 ; 

777(771774) (3 1 1 1 2] ( (2 1 1 1 2] + 2777^) 777 (2771 ) (27/4) (2 1 1 1 2] (3 1 1 1 2] 



F, 
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(2|1|3]S23 (23)(2|1|3]^ 



+ 
+ 



777(3774) (771 1 1|2](3|1|2] 777(377i)(37/4)(3|l|2](27n=^-(2|l|2]) 2777=^ (2771) (3774) (3|1|2] 

2(2|l|2]si2 ^ 2(23)(2|1|2]2 + (23)(2|1|2](2|1|3] 

777(27/i)(3774)(3|l|2]=^ 777(27;i)(27/4)(3|l|2]^ (2(2|1|2]^ + 2(2|1|3](3|1|2] + (2|1|2]523) \ 
(23)(2|1|2]2 2(23)(2|1|2]2(2|1|3]S23 J 

777(37;4)(3|1|2] ((277i)(3|l|2] - (37/i)(2|l|2]) 2777(77i774)(2|l|2](3|l|2] 777(771774) (3|1|2]^ 

2(23)(2|l|2]si2 (2|l|3]s23 2(2|l|2]s23 

777(277i)(3|l|2]((2774)(3|l|2]- 4(3774)(2|l|2]) 777(2771) (2774) (3|1|2]^ 777(377i)(3774)(3|l|2] 
2(23)(2|1|2](2|1|3] (23)(2|l|3]s23^2 + (23)^23^2 

2777(277i)(27/4)(3|l|2]((2|l|2] + 27772) 

(23)(2|l|3]2/32 17 

(94) 



^(r/in(r/44VP(l^2+,3-,4) = 

Af^ (1+ ,2+,3~,4^^) X (^(s23-2m^)-^^l|23|4 (2|l|2]/32 -^3?2|3|41 + ^2|1|2]^2;12 " ^h;m 

+ (-((2|1|2] + 2m^)[23]/5j^-J, + ^-^^^^^I^,,,,,, - m^[23]I-,,,,, + ^^^|^i^2;i2 
_(SM(j,m ,A m(277i)(3|l|2]((2774)(2|l|2] + 2(3774)(2|l|3]) e-2e.. 

(23) ^ 2;23 ^J^ X (2|1|2](2|1|3]2 ^4;1|2|3|4X 

/ m^(2r74)(3ryi)(3|l|2]2 m(3|l|2][23](3r7i)(3774)((2|l|2] + 2m^) m=^(3r?i)(3r?4)(3|l|2][23] 

^ V (23)(2|l|2]2/32 + (2|1|2]2 ^ 2(2|1|2]2 

m(277i)(3r^4)(3|l|2](2(2|l|2]2 + g23((2|l|2] + 2m^)) m(27yi)(2r74)(3|l|2]2 ((2|1|2] + 2m^) 

(23)(2|1|2](2|1|3] + (23)(2|l|3]s23/32 

m(2r;i)(3|l|2]2((2r74)(3|l|2] + (3r;4)(2|l|2]) 3m(2r;i)(3r74)(3|l|2]2 ((2|1|2]2 + 51^^23) 

(23)(2|l|2]s23/32 + (23)(2|l|2]2s23/32 

2m(3r/i)(3^74)(3|l|2]((2|l|2] + 2m") m3(3r]i)(3r74)(3|l|2] (2(2|1|2] + S23) 



(23)(2|l|2]/32 2(23)(2|l|2]2/32 
m{vm){msi2{{2\l\2f - s,2S23) \ ^ m(3r/4)(3|l|2] (2(2r/i)(3|l|2] + Sm^jSy,)) 

(2|1|3](2|1|2]2 ^+^4;l|2|3|4 ^33) 



/5/,r'^"(3//i)(3//i)(3|l|2] 2rrr^(//i//0 (3|l|2],si2 m3(3//i)(3//i)(3|l|2] (8( 



+ 3;i2|3|4\^ (23)(2|1|2] (2|1|2](2|1|3] (23)(2|l|2]s23/32 

2m^(77i774)(3|l|2](((2|l|2] + 2m^)(2|l|2] + (2|l|3](3|l|2]) 2m3(2r;i)(2r/4) (3|1|2]2 

(2|l|2](2|l|3]s23/32 ^ (23)(2|l|3].23/92 

2m3(27yi)(3r/4)(3|l|2] (2(2|1|2] + ^23) 4m(277i)(3r/4)(3|l|2]2 ((2|1|2] + 2m^) \ 
(23)(2|l|3]s23/32 (23)(2|l|2]s23/32 ) 

,, m(3|l|2]2((2r^4)(3r7i) + (2r;i)(3r74)) /m(2r/i)(3r/4)(3|l|2][23] 

V-^2;23 ^) /OqXea /Q2 -'3;2|3|41 I 



■''23 J 



(23).23/32 (2|i|3] 
2(^1^4) (3|l|2]2m3 2m3(r/ir/4)(3|l|2] ((2|1|2] + 2m2) 6771^(3771) (3774) (3|1|2] [23] 

(2|l|2]s23/32 ^ (2|l|3]s23/32 (2|l|2]2s23/32 

r7i(277i)(3774)(3|l|2] (2r7i2(2|l|2]2 - (2|1|3](3|1|2]((2|1|2] + 4771^)) 



(23)(2|1|2]2(2|1|3]^2 

77i(277i)(3|l|2] ((2774)(3|l|2]5i2 - (774|1|2](23)(2|1|2]) \ / 2771(771774) (3|1| 2] 512 

(23)(2|l|2](2|l|3r )^ (2|1|3]S23 

777^(377i)(3774)(3|l|2] _ 77i(277i)(2774)(3|l|2]2 _ 777(277i)(3774)(3|l|2]2 777(2771) (3774) (3|1|2] \ 
2(23)(2|l|2]si2 2(23)(2|1|2](2|1|3] 2(23)(2|l|2]si2 (23)(2|1|3] ) 

, ^, ^ 777(77i774)(2|l|2](3|l|2] , 777(77i774)(3|l|2]((2|l|2] + 2777^) 

^'^^^ (2|1|3]S23 + (2|l|3].23/3^ 

777(3|1|2] (2(2|1|2] + 3523) ((377i)(3774)(2|l|3] - (277i)(2774)(3|l|2]) 

2(23)(2|l|3]s23/32 

777(277i)(3774)(3|l|2] (2(2|1|2] + 523) \ 77i(3|l|2]^((77i774)(2|l|3] + (277i)(2774)[23]) 
(23)(2|l|3]s23;52 ; 2(2|l|2](2|l|3]s23 

r»(3/A)(3|l|2]((3//i)(2|l|2]-(2//i)(3|l|2]) 
2(23)(2|l|2]si2 

(95) 
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The primitive amplitudes amplitudes contributing to A4.3 contain bubble, tadpole and 
rational terms that cancel when forming the complete amplitude. We do not list these terms 
explicitly in the following expressions. 



-t{Vil'){v^^')AP (l+,2+,4t 3+) = 

'm^rim) ((3|1|3]2 - (2|l|2]s23) m3(3r/i)(3r74)[23](2|l|2](3|l|3]\ 



+ -^4;l|2|4|3 X 



(23)2 2(23)2(3|1|2] J 

(2(r;ir/4)(3|l|2] + {3r],){3r]4.)[23]) 



+ ((2|1|2] 73;12|3|4 + (3|1|3] 73;13|2|4) X 2(23)2(3| 1 12] 

,(/n\,\num ./qiiiqU- ^ m'^ (2far?4) (3 1 1 12] - (3r?,) (3r?4) [23]) 

+ l\^|-LMJ^3;12|3|4+ W|-L|>iJ^3;13|2|4j X 2(23)2(3| l|2] 

,(rm _jm ^. ^^[23] ((277i) (3^74) + (377i) (27^4)) 
+ l-'3;12|3|4 -'3;13|2|4j X ^^3^2 

/^3Tl2|3|4 , ^37l3|2|4\ m3(277i) (2^4) [23] (3| 1 12] 

1,(2|1|2] ^ (31113]; (23)2 
+ bubbles, tadpoles and rational terms (96) 



-i{mi'){iiA^')A^i:^ (1+2+, 4, 3-) = 

^m((r;ir;4)(3|l|2] - (3r7i)(3r74)[23]) (-(2|1|3](3|1|2] + m2s23) 

+ -'4;l|2|4|3 X 
+ 



2(2|1|3] 

m^(2|l|2](3|l|3] ((2^7i)(3774)[23] + (7717/4)^12) 3m3(3r/i)(3^74)(2|l|2](3|l|3] 

(2|1|3]2 2(23)(2|1|3] 
2m=^(277i)(3774)(2|l|2]2(3|l|3] ^ m3(277i)(2774)(2|l|2]2(3|l|3]2 



(23)(2|1|3]2 2(23)(2|1|3] 

+ (-^3^12|3|4 ~ -^3;12|3|4) X ^ 



m=^(77i774)(2]l]2]si2 ^ m3(277i)(2774)(2]l|2]2(3]l|3] 



+ 



(2|1|3]2 2(23)(2|1|3]3 
m(2|ll2] ((771774) (3|1|2] - (3771) (3774) [23]) 3m=^(377i)(3774)(2|l|2] 

2(2|1|3] 2(23)(2|1|3] 
m3(277i)(37/4)(2|l|2][23] 27^3(2771) (3774) (2|1|2](3|1|3]^^ 



(2]1|3]2 (23)(2|1|3]2 ; 

+ (-^3?13|2|4 - -^3;13|2|4) X ^ 



7n3(r7i774)(3|l|3]si2 ^ 7n3(277i)(2774)(2|l|2](3|l|3]^ 



(2|1|3]2 2(23)(2|1|3]3 
7n(3|l|3] ((771774) (3|1|2] - (377i)(37/4)[23]) 37n-3(377i)(3774)(3|l|3] 

2(2|1|3] 2(23)(2|1|3] 
7n^(277i)(3774)(3ll|3][23] 27n^(277i)(3774)(2|l|2](3|l|3] \ 
(2|1|3]2 (23)(211]3]2 ; 

+ bubbles, tadpoles and rational terms (97) 



+ 
+ 
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C. Primitive amplitudes for qq tt 



i{r^^l'){r^,4')Al^\li,2l3-,4t) = 
if (1^2+ 3-,4±) X {^^^^/4;i|2|3|4-f Wi + (2|^ 

3m2(2(2|l|2] + S23)\ , 8^ 23 12m2(2(2|l|2] + S23) 



-{S23+ j /3;l|23|4 - 'h^m + y + ^^^^4 

7 (2(2|l|2] + g23) 6m^(2(2|l|2] + g23) \ f 
G 2523/52 + sl,(3^ J '-''^ 

m((2r^i)(2r/4)[23] + far?4)(2|l|3]) / ,..3, ,.12, 

2(2|1|3]2 |\2|1|2J -/4;1|2|3|4- \2|1|2J i3;2|3|41 

- ^^^^^/3;12|3|4 + (S12S23 - 2m2(2|l|2]) /3;i|23|4 
•523 

2 ,f o^^ /6m3((377i)(3?74)(2|l|2]2 + (2^i)(2?74)(3|l|2]2) 



+ (^^3;l|23|4 + -^{h;23 " 2)^ X 



+ 



6m^(2r^i)(37y4)(3|l|2](2(2|l|2] + ^23) _ 3m=^(3r/i)(3r/4) \ (98) 
(23)523/34 2(23)/34 

/ m(2r;i)(2r74)(3|l|2]((2|l|2] + 2m2) 2m(2^i)(3774)(3|l|2] 
V (23)(2|l|3]s23/52 ^ (23).23/32 

m(3r7i)(3r74)gi2 \ / m ((277i)(2774)(3|l|2] + (3r7i)(3774)(2|l|3]) [23] 

(23)s23/32 y' '='I''I'V 2(2|1|3] 
m3 (3771) (3774) (4(2|1|2] + S23) m^2r]i){2r]4){3\l\2] (2(2|1|2] + S23) 



^-^2;23 ~ 2^ 



2 



2(23)523/52 (23)(2|l|3]s23/5' 
4m=^(277i)(3r^4)(3|l|2]^ , ^ /m((377i)(3r^4)(2|l|2]2 + (277i)(2r^4)(3|l|2]2) 

-'2:23 I 



(23). 23/32 ; '^^'v (23)^23^2 

m(27?i)(3r74)(3|l|2](2(2|l|2] + S23) m^3vi){3v4)\ ^ ^ f m{27i^) {2r],) {3\1\2]' 

+ -^2:12 I 



(23)4,/32 (23)s23/32 ; '^^'V S23(23)(2|l|2] 

m(2r/i)(2r74)(2|l|2](3|l|2] 771(37^4) ((2r/i) (3| 1 12] + m2(377i)) \ 
S23(23)(2|l|3] (23)(2|1|2] J 

-z(7;il^)(r?44Vr](l+ 2+ 3-,4±) = 

- if (1+ 2+ 37, 4±) X (^S23/3;2|3|41 + (^23 - 277l2) 123,4 

3^ / 2(2|l|2] + g23 \ , , , 2(2|1|2] + S23 

+ ;7-'2;23 + J- H /Q2 "^2:23 + ^2;m " 4 



2^'^^ V 2S23/32 ; ^'^^ S23/32 

(3771) (3774) , 7n(27yi)(37/4)(3|l|2](2(2|l|2] + S23) 



/™ ^^ ^ f m^3rii){3ri4) 
- (^2;23 - 2) X ^23).23/32 + 



(23)523/32 

m((377i)(3774)(2|l|2]2 + (277i)(2774)(3|l|2]2 )' 
(23)^23^2 
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(99) 



D. Fermion loop amplitudes 



Below we list the heavy fermion loop corrections for an arbitrary mass m^/. The light 
fermion loop contributions can be obtained by taking mn 0. 

-z(r/il^)MVf (l,+ ,2+,3+,4±) = 

2m(fary4)(2|l|2]-(2ry,)(3ry4)[23]) / F^" + \ ] (100) 

(23)3[23] V "^^^^3;2|3|41 + ^^H^2;23 + g'523 I , liUUj 

-^Af (1+ 2+,3-,4) =0, (101) 
-^Am (1+ 2+ 3-,4±) = 

-|4°^(1+ 2+ 3-,4±) ((2^ + 1) F-%+I---l^ . (102) 



IX. NUMERICAL RESULTS 

The amplitudes in the previous section have been implemented into an efficient library 
for evaluation of the colour and helicity summed interference with tree level. In order to 
minimise the number of independent spinor products a specific choice of ?7i = 774 = P2 was 
made. This raises a number of issues since the symmetries between the helicity ampli- 
tudes are broken, however it is straightforward to generate all the necessary configurations 
automatically from those presented here. 

For illustrative purposes we present numerical values for the unrenormalised amplitudes 
in the FDH scheme with the strong coupling set to one. We choose a generic phase-space 
point for the momenta as follows: 

Pi = (-Vs + m2,yi,0,o) , (103) 

P2 = V s + m'^ (1, sin 6, cos 6 cos 0, cos 6 sin 0) , (104) 
P2 = V s + m? (1, — sin 9, — cos 9 cos 0, — cos 9 sin cp) , (105) 
Pi = (^-Vs + m2,-yi,0,o) . (106) 

(107) 

The mass of the heavy quark is m = 1.75. Numerical results at the point s = 1, 6* = |, 
(f) = J, Nf = 5, A'^^^ = 1, A^c = 3 are given in Table [11 The interference with the tree 
level amplitude summed over helicity and colour is given in Table [III at two values of the 
renormalisation scale, /i^ = and = Am^. 

Our results have been cross-checked against previous calculations in the literature. At 
the amplitude level we find full agreement with numerical results of 19 1. We have also 



checked against the analytic results of ref. 16|] up to 0(e°). We find full agreement with 



the implementation of the results of ref. j46|] into MCFM [47]- 

The FORTRAN program used to generate the results of Table [III evaluates the interfer- 
ence of one-loop amplitudes with the tree-level, summed over helicity and colour, in 43/is 



A description of the reference vectors used to verify the resuhs of ref. is given in Appendix [Bl 
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Primitive Amplitude 




(l+,2+,3+,4±) 


0.055220794 + 0.014807839i 


2+,3-,4±) 

4 V t ' ' ^ t 


0.062949503 + 0.14075936i 


aW(1+ 2+ S+ 4+^ 


50453481 -1- 3385402? 


Ai^\(l+ 2+ 3+ 4+1 


021 1 00789 — 1 2891 5fi37 


^[,^'(1,+ 2+ 3+ 4^) 


—0 01309239 + 028932428z 


Af\(l+ 2+ 4+1 


—0 00391 70375 — 04534929? 


yl..nfl+ 2+ 3+ 4^1 


63758829 - 1 1392369? 


^i^Vl,+ ,2+,3-,4+) 


1.1550236 + 1.3088169i 


Af^(l+,2+,3-,4±) 


2.5565516 + 1.0995254i 


4''kl^2+,3-,4±) 





4-'''(l+,2+,3-,4±) 





A4;3(1^2+,3-,4+) 


11.18323 - 1.475571i 


<\l+,2± 3-,4+) 


0.85072714 + 0.25682619i 


2±,3-,4±) 


8.0971525 + 3.2796876i 


4^'^J(l+,2+,3-,4+) 


15.65914 - 4.9530347i 


<'(1^2± 3-,4+) 


1.3699007 + 1.7416922i 


Alt\(-[+ 0+ Q- 4.+ \ 


0.24709861 + 2.0187408i 



TABLE I. Numerical values for the individual primitive amplitudes. The reference vectors are 
chosen as ??i = (3, 2, 2, 1) and 7/4 = (3, 2, 1, 2), /i/j = 2m. 



Channel 


1 


1 

e 


6° 


99 ifiR = m) 


-882.7183832 


1406.029038 - 1915.339983i 


2811.35321 + 1478.791625i 


99 tt {^R = 2m) 


-882.7183832 


182.3215209 - 1915.339983i 


3912.313922 - 1176.433394ii 


qq tt {nR = m) 


-82.70769231 


313.1028325 + 73.26833668i 


152.9616128 - 326.7068046ii 


qq -> ti (nR = 2m) 


-82.70769231 


198.4456251 + 73.26833668i 


507.5399839 - 225.1353226ii 



TABLE II. Numerical values for the interference between virtual and tree level amplitudes summed 
over helicity and colour, 2'^^j^A^I\a^^'^]*. 



for the gluon fusion channel and 13/is for the quark annihilation channel. Roughly 60% 
of this time is spent on the evaluation of the scalar integrals. The code was compiled us- 
ing gf ortran with optimisation level -02 and evaluated on a 2.93 GHZ Intel Core 13 530 
CPU. The FORTRAN code BSYpptt and the FORM files used to produce it are available 
from http : //www . nbia . dk/badger . html . Both colour and helicity summed results and the 
individual primitive amplitudes are included. 



X. CONCLUSIONS 

In this paper we have used the newly developed techniques of generalised unitarity to 
compute compact analytic representations of all helicity amplitudes relevant for heavy quark 
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production at hadron colliders. Compact tree-level input was generated via BCFW recursion 
relations and the coefficients of the scalar integral computed via a purely algebraic procedure. 
A fully automated Feynman diagram approach was used to produce compact forms for the 
tadpole and rational terms. 

The calculation was performed in the spinor-helicity formalism with a completely general 
representation for the heavy quark spinors. The final amplitudes are expressed in terms of a 
relatively small set of spinor products. The analytic forms of the helicity amplitudes allow us 
to take a new look at the structure of the one-loop helicity amplitudes. The all-plus helicity 
configuration in the gluon channel takes a remarkably simple structure. In comparison to 
the well known MHV structure in massless amplitudes it is expected that such a simplicity 
persists at higher multiplicity. The most notable feature of this new representation is the 
cancellation of bubbles and rational terms in the sub-leading colour contribution to gg — )■ tt 
which is expected to be related to the stronger UV constraints on this sector. Together 
with similar features found in other analytic computations |^], this may motivate future 
investigations. 

The final results yield a flexible implementation suitable for computations of spin corre- 
lations and decays in the narrow width approximation. 

We have demonstrated that such techniques provide a feasible method to calculate an- 
alytic one-loop amplitudes with full mass dependence and serve as a solid base for future 
implementations of higher multiplicity amplitudes. 



ACKNOWLEDGMENTS 

We are particularly grateful to Keith Ellis, John Campbell, Peter Uwer and Sven-Olaf 
Moch for helpful discussions and feedback throughout this project. The work of VY was 
supported by DEC SEB-TR-9 and by LHCPhenoNet network PITN-GA-2010-264564. The 
work of SB has been supported in part by Danish Natural Science Research Council grant 10- 
084954. The work of RS was supported in part by grant GK 1504 of the DEC Graduiertenkol- 
leg. 



Appendix A: Tree-level amplitudes 

Eor completeness we present the tree-level amplitudes relevant for our computation using 
the spinor-helicity formalism described in section IIIII 



1. On-shell Three-point Vertices 

The independent on-shell three-point vertices are: 



^^3(1^,2+, 3^ 



^A3(l+,2+,3^ 



+ 



(e|l|2] mfar^a) 

(e|l|2](1^3^) 
(^2) m ' 



(A2) 



(Al) 
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-i^3(l^.2^3^) = -W^. (A3) 

" ' ' {{2> (li'ifa) ^ ' 

2. gi^i — )• Tree Amplitudes 

For the gg — ?■ it channel with adjacent fermions we obtain 

lA, [i, ,Z ,6 ,4,j - ^23)(2|l|2](r/il^)(r/44^) ' ^ ^ 



.^(0).^+ 2+ 3- 4+) = "^(3|l|2]((r/ir^4)(3|l|2]-[23](3r/i)(3r^4)) . 
' ^ *' ' ' *^ S23(2|l|2](r7il^)(r744^) ' ^ ^ 



For the sub-leading colour contributions we also make use of compact forms for the case of 
non- adjacent fermions: 

_,,4(o)n+ 2+ 4+ 3+) - ^'fa^4)[23]^ 

' ^ ' ^" (r/il^)(r/44^)(2|l|2](3|l|3]' ^""'^ 

2+ 4+ 3-) = "^(3|l|2](fa^4)(3|l|2]-[23](3r/i)(3r^4)) 
' ^ *' ' ' (r/il^)(r/44^)(2|l|2](3|l|3] ' ^ ^ 

The other fermion helicity states can be obtained via the relation given in (123|) . 



as, 



3. qq — > tt Tree Amplitudes 

There is only one independent helicity amplitude in this channel which can be written 

m((r7i3)(r74|4|2] + (r/43)(r7i|l|2]) 



-iA[ 2^, 3^ , 4^ 



(?7il^)(r744^)s23 



(A9) 



Appendix B: Conversion to Four Component Dirac Spinors 



We note that the formalism of equations f[T9l) and fl20l) can be connected with a more 
conventional approach to massive solutions of the Dirac equation by choosing a specific 
reference frame. For a massive four-vector = {E,Qi,Q2,Q3) with = m? we first 
define: 



Q+ = E + Q3, 
= Qi + iQ2, 



Making a choice of 



1 



2(g_ + m)2 



Q^=E-Q,, 
Q^ = Qi- 1Q2. 



( + mf + Ql + Ql\ 
2Qi(g_ +m) 
2Q2{Q- + m) 
\-{Q.+mf + Ql + Qll 



(Bl) 
(B2) 



(B3) 
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then yields the following four dimensional representations: 





VE+m 
\ / 



VE+m 



VE+m 

y/E + m 




( \ 

\/E^m 



VE+m 

Qz 



VE+m 



VE+m 



VE+m 





(B4) 



(B5) 



\VE + m) 

This allows a simple way to compare analytic results with numerical ones in the literature. 
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